Conservative 3+1 General Relativistic Boltzmann Equation 
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We present a new derivation of the conservative form of the general relativistic Boltzmann equa- 
tion and specialize it to the 3+1 metric. The resulting transport equation is intended for use in 
simulations involving numerical relativity, particularly in the absence of spherical symmetry. The 
independent variables are lab frame coordinate basis spacetime position components and comoving 
frame curvilinear momentum space coordinates. With an eye towards astrophysical applications — 
such as core-collapse supernovae and compact object mergers — in which the fluid includes nuclei 
and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we 
examine the relationship between lepton number and four-momentum exchange between neutrinos 
and the fluid. 

PACS numbers: 95.30.Jx, 05.20.Dd, 47.70.-n, 97.60.Bw 



INTRODUCTION 



Neutrino transport is a necessary ingredient of core- 
collapse supernova simulations [1-7]. Determining the 
fate of the stellar material — for instance, does an explo- 
sion happen, and if so, how? — requires calculation of the 
four-momentum and lepton number exchange between 
the fluid (which includes nuclei and / or nuclear matter at 
finite temperature) and the neutrinos that stream from 
and through it. For the purpose of studying the explosion 
mechanism, we take the traditional approach and con- 
sider only massless neutrinos described by classical dis- 
tribution functions (phase space densities) /(t, x, p) 8]. 
By way of introduction to the Boltzmann equation satis- 
fied by f(t, x, p), we discuss at some length the choices of 
reference frames for the spacetime and momentum space 
coordinates; and more briefly, conservative formulations 
and efforts to maintain consistency of lepton number and 
four-momentum exchange with the fluid. 

The form of the Boltzmann equation governing dis- 
tribution functions f(t, x, p) depends on the choices of 
spacetime and momentum space coordinates. For an 
overview of possible effects, consider the Boltzmann 
equation expressed as 



df 

at 



vV x / + F-V p / 



9/ 

at 



(i) 



collisions 



in which the number density of particles at position x 
with momentum p changes in time because (a) they 
stream away from x with velocity v = dx/dt = p/E(p) 
(where E(p) is the particle energy); (b) their momenta 
are gradually altered by long-range (i.e. electromag- 
netic and/or gravitational) forces F = dp/dt; and (c) 
their momenta are abruptly altered by point-like colli- 
sions. Naturally the precise expression of the gradients 
depends on the choice of coordinate systems. But less 
trivially, curvilinear spatial coordinates and rotating ref- 
erence frames respectively would induce 'ficticious' and 



Coriolis force contributions to F. And as discussed fur- 
ther below, one may wish to reckon the components of x 
and p with respect to two different reference frames in rel- 
ative motion; this further complicates the left-hand side 
with momentum- (component)-changing contributions to 
F, associated for instance with angular aberrations and 
Doppler shifts. Sometimes called 'observer corrections,' 
these terms involve the relative velocity of the two frames, 
and — in the case of accelerated relative motion — its gra- 
dients. On the importance of such effects, see for instance 
Ref. 0. 

Lindquist [l(| introduced the use of distinct coordinate 
systems for spacetime and momentum space in the con- 
text of general relativity, through use of an orthonormal 
basis or 'tetrad.' Given a spacetime metric in a coordi- 
nate basis (also called a 'natural' or 'holonomic' basis), an 
orthonormal basis can be chosen at each spacetime point, 
in terms of which the metric locally takes the Minkowski 
form. Lindquist began by writing the general relativistic 
Boltzmann equation in terms of such an orthonormal ba- 
sis. He then re-expressed all time and space derivatives — 
df/dt and V x / in Eq. ([T]), as well as the metric and 
tetrad derivatives entering F — in terms of the original 
coordinate basis. However, all momentum components, 
and the momentum space gradient V p / in Eq. (JXJ) , were 
left in terms of the orthonormal basis. Thus the indepen- 
dent variables were the spacetime position components in 
the coordinate basis and the momentum components in 
the orthonormal basis. 

What purposes did this use of distinct coordinate sys- 
tems for spacetime and momentum space serve? An ob- 
vious advantage of an orthonormal basis in momentum 
space — not mentioned by Lindquist, and rarely (if ever) 
noted explicitly, perhaps because of its obviousness — is 
that, as a local Minkowski frame, it facilitates direct 
importation of flat-space results for particle interaction 
rates into the collision integral. As for Lindquist's pur- 
poses, one was to make contact with familiar expressions 
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from classical transport theory, and in some respects the 
use of an orthonormal momentum basis helped with this. 
He also noted that the 'tangent bundle' (spacetime plus 
the tangent space at every point of spacetime) is a natu- 
ral way to conceptualize phase space in general relativity; 
in this context, in which momentum space is identified 
with the (flat) tangent space at each spacetime point, an 
orthonormal basis for momentum space also seems nat- 
ural. But for Lindquist, a primary benefit of the 'tetrad 
formalism' was not necessarily the orthonormality per 
se, but the mere fact that introduction of a second co- 
ordinate system allowed simultaneous exploitation of the 
differing symmetries of the spacetime and the associated 
momentum space he considered (i.e. spherical symmetry 
in position space, vs. axial symmetry about the radial 
direction in momentum space). 

While Lindquist used different bases for position and 
momentum space, he did not mention the possibility that 
these coordinate systems could be in relative motion. 
Having in mind spherically symmetric gravitational col- 
lapse, it is perhaps not surprising that in specializing his 
Boltzmann equation he considered only comoving coor- 
dinates in both position space and momentum space. 

The possibility of relative motion between the position 
space and momentum space coordinate systems appar- 
ently arose first in the context of special relativity. Cas- 
tor [ll| specialized Lindquist 's comoving frame results in 
spherical symmetry 10] to flat spacetime, keeping only 
terms valid to 0(V/c) in the fluid speed V. Regard- 
ing this, Mihalas [12|, LL3| noted that the collision inte- 
grals were most easily evaluated in the comoving frame, 
for this is where the distributions of particle species in 
equilibrium (e.g. in the fluid, but also the transported 
species at high optical depth) were isotropic. But he ex- 
pressed concern that it is only in an inertial frame that 
"synchronism of clocks can be effected," and that use 
of inertial-frame time and space coordinates "obviates 
the need to develop a metric for accelerated fluid frames 
(Castor 1972) which in general can be done only approx- 
imately." If anywhere, such concerns ultimately could be 
relevant only within the confines of special relativity, but 
they nevertheless motivated Mihalas to "leave both the 
space and time variables in the inertial frame" (emphasis 
in original; while there are an infinite number of inertial 
frames, "the" inertial frame he used was a 'lab frame' at 
rest relative to the origin of some spherically symmetric 
physical system). In Ref. [l2j, consistent with a focus 
on special relativity, he eschewed Lindquist's general rel- 
ativistic machinery. Instead, after writing down the lab 
frame transfer equation, he used the Lorentz transfor- 
mation properties of the radiation variables to obtain a 
radiation transport equation in spherical symmetry, valid 
for < V/c < 1, whose independent variables were the 
lab frame time and radius and the comoving frame par- 
ticle energy and direction angle. 

The most important motivation to reckon spacetime 
coordinates in an Eulerian frame was one beyond the 
focus of Ref. [l2j on spherical symmetry: in multiple 



spatial dimensions, mesh-based treatments of the fluid 
with which the particles interact are often formulated in 
an Eulerian frame, because such phenomena as differen- 
tial rotation, convection, and turbulence result in mesh 
stretching and entanglement in a Lagrangian approach. 
Using to some extent general relativistic perspectives, flat 
spacetime transport equations in full dimensionality, ex- 
pressed in terms of lab frame spacetime coordinates and 
comoving frame momentum space coordinates, were ob- 
tained by Buchler [3, [l5| and by Kaneko et al. [HI to 
Q(V/c), and by Munier and Weaver for < V/c < 1 

EE Eg. 

Returning to Lindquist [Io| as a point of departure, a 
fully general relativistic Boltzmann equation, expressed 
in terms of lab frame coordinate basis spacetime coordi- 
nates and orthonormal comoving frame momentum vari- 
ables, was given independently by Riffert [l9[ and by 
Mezzacappa and Matzner [2(|. While they did not ex- 
press it quite this way, the basic point in connection with 
Lindquist's approach is this: the transformation from a 
coordinate basis to an orthonormal basis — that is, the lo- 
cal diagonalization of the metric to the special relativistic 
Minkowski form at each point in spacetime — is unique 
only up to rotations in spacetime, which includes boosts 
as well as spatial rotations. In other words, even if one's 
coordinate basis is some sort of Eulerian lab frame, one is 
still free to choose Lindquist's tetrad of orthonormal basis 
vectors to be one carried by a Lagrangian observer, one 
comoving with the fluid and therefore (in general) under- 
going accelerated motion relative to the lab frame. This 
can be conceptualized as two sequential transformations, 
first from the lab frame coordinate basis to an orthonor- 
mal lab frame basis, followed by a Lorentz boost [20|; or 
directly as a single transformation, in which the time- 
like orthonormal basis vector is the fluid four-velocity 
[l9j (see also Ref. 21 1 and, in connection with comoving 
spacetime coordinates, Ref. [22]). Mathematically, accel- 
eration of the Lagrangian observer appears as velocity 
derivatives in the comoving frame connection coefficients 
entering F in Eq. ([1]). 

The presentation of Mezzacappa and Matzner [2(| in 
particular, couched as it is in terms of the geometric con- 
cepts of the 3+1 approach — and conceived with a view 
towards multidimensional simulations in which spacetime 
evolution, hydrodynamics, and Boltzmann neutrino radi- 
ation transport are all treated in full general relativity — 
can be regarded as the inspiration for the extension given 
in the present work. Full elaboration of the general rel- 
ativistic Boltzmann equation in lab frame spacetime co- 
ordinates and comoving frame momentum variables was 
only given in spherical symmetry in Refs. [l9l [2fjj . and 
the complexity was already nontrivial in that case. Be- 
yond the mere multiplication of dimensions, contribut- 
ing to the complexity in the multidimensional case is a 
non-diagonal spatial metric that effectively must be diag- 
onalized in the transformation to an orthonormal frame. 
In the face of an inherently complicated equation, we 
strive to minimize pain and maximize utility to readers, 
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by aiming for an appropriate balance between maintain- 
ing simplicity and clarity while keeping enough detail to 
make numerical implementation straightforward. 

We now briefly discuss conservative formulations and 
efforts to maintain numerical consistency of lepton num- 
ber and four-momentum exchange with the fluid. In a 
conservative formulation, the left-hand side of Eq. ([T]) 
is rewritten as a phase space divergence — that is, as 
the sum of spacetime and momentum space diver- 
gences. Conservative formulations aid efforts to numeri- 
cally maintain fidelity to global conservation laws (such 
as lepton number and energy) at modest resolution by 
transparently relating volume integrals to surface inte- 
grals. A conservative formulation of the Boltzmann equa- 
tion in spherical symmetry to 0(V/c) was given by Mez- 
zacappa and Bruenn (23|. Among other important dis- 
cretization considerations, they generalized to the Boltz- 
mann case a scheme for an explicit update of / from 
the energy derivative term that was consistent with both 
lepton number and energy conservation [24| (originally 
developed by Bruenn [251 ] for the case of flux-limited 
diffusion). Still in spherical symmetry but in full gen- 
eral relativity, Liebendorfer et al. [26[ generalized this 
scheme to the angle derivative terms; and in order to 
avoid operator splitting of the observer corrections, con- 
verted these energy and angle updates into finite differ- 
ence representations suitable for inclusion in the implicit 
solve used for the other terms in the Boltzmann equa- 
tion. Liebendorfer et al. [26| also deduced discretizations 
of various terms in the number-conservative Boltzmann 
equation that provided consistency with the lab frame 
energy conservation law. Cardall and Mezzacappa (27| 
gave number- and four-momentum-conservative formula- 
tions of the Boltzmann equation in full general relativity 
and in the full dimensionality of phase space, but not 
in a 3+1 form useful in numerical relativity. First tests 
of a full 3D+3D Boltzmann solver have been reported 
[1, HI| i but without all the terms that appear as a result 
of fluid motion, which provide challenges to simultaneous 
energy-momentum and lepton number conservation. 

Due to the reduced computational burdens incident to 
integrating out the momentum space angular variables, 
moments approaches have received recent attention as 
the field turns towards simulations that are 3D in posi- 
tion space. In a conformally flat ray- by-ray approxima- 
tion, Miiller et al. [29[ developed a scheme that renders 
non-conservative (in a sense discussed in Ref. (30[) energy 
and momentum equations consistent with a conservative 
number equation, and have deployed it in 3D position 
space simulations retaining energy dependence [311 ] . (An- 
other energy-dependent ray-by-ray simulation in 3D posi- 
tion space has been performed with a different approach, 
the 'isotropic diffusion source approximation' [321] .) Shi- 
bata et al. [33[ presented conservative energy and mo- 
mentum angular moment transport equations in the 3+1 
formulation of general relativity. These have been imple- 
mented in simplified, energy-integrated ('grey') form in 
the 3D position space simulations of Kuroda et al. f34j . 



In order to more fully address simultaneous energy and 
lepton number conservation in moments formalisms, the 
energy derivative term in the Shibata et al. [33| equa- 
tions was more fully elaborated into a useful 3+1 form by 
Cardall et al. [35j , who also considered in detail the term- 
by-term cancellations that must occur for consistency be- 
tween conservative four-momentum transport equations 
and the conservative number transport equation; see also 
the partially relativistic limits discussed by Endeve et al. 
30]. 

In this paper we build on these developments by pre- 
senting conservative formulations of the Boltzmann equa- 
tion in the 3+1 approach to general relativity. We 
present a new derivation of the number-conservative for- 
mulation in Sec. [TTJ The independent variables are lab 
frame coordinate basis spacetime position components 
and comoving frame curvilinear momentum space coor- 
dinates. In Sec. IIIII we use 'Eulerian decompositions' of 
the tetrad and its inverse in specializing to the 3+1 met- 
ric, and examine the relationship between the number- 
conservative and four-momentum conservative reformu- 
lations. We collect and discuss our results in Sec. IIV1 
including overview tables of the many variables appear- 
ing in the formalism. 

We conclude this introductory section by specifying 
some conventions we employ. We take h = c = 1 and use 

the metric signature h H — h Greek indices run from 

to 3 and latin indices run from 1 to 3. We denote 
lab frame coordinate basis components with unadorned 
indices; components measured in the orthonormal co- 
moving frame with hatted indices; and momentum space 
components in a curvilinear momentum space coordinate 
system with indices adorned with a tilde. 

II. CONSERVATIVE FORMULATION 

After specifying phase space coordinates, we present 
a new derivation of the conservative formulation of the 
general relativistic Boltzmann equation. See Ref. [27| for 
an alternative derivation using exterior calculus. 

A. Phase space coordinates 

The general relativistic Boltzmann equation [ToL 
I20I I3H [37} equates the directional derivative of the dis- 
tribution function along the phase flow to a collision in- 
tegral: 

That is, the change in / along a phase space trajectory 
with afflnc parameter A is equal to the phase space den- 
sity C[f] of point-like collisions that add or remove par- 
ticles from the trajectory. The phase space measure is 
defined in such a way that / and C[f] are both invariant 
scalars. 
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For practical computations it is necessary to introduce 
phase space coordinates, which we initially take to be 
spacetime position coordinates x p and momentum space 
coordinates p % measured in an orthonormal frame comov- 
ing with the fluid. We denote components in the comov- 
ing frame with hatted indices. We take only the spatial 
momentum components p l as independent variables, due 
to the mass shell relation 

m 2 = - PA pA = -ViXuPh" , (3) 

where m is the particle mass and (rjfip) — diag[— 1, 1, 1, 1] 
is the Minkowski metric. In terms of these phase space 
coordinates, Eq. (JTJ becomes 



dx^ 1 df dp 1 df 



dX dx^ dX dp' 



I = C[f]. 



(4) 



For electrically neutral particles, the geodesic equations 
specifying particle trajectories are 



dx» 
~dX 
dp» 
~dX 



P" 



so that 



(5) 
(6) 

(7) 



(For electrically charged particles the coupling to the 
electromagnetic field can come along for the ride in all 
that follows; see e.g. Ref. [27}.) As we shall see (cf. 
Eqs. [[I0 |) -(|T2 ]l below), when spacetime is curved, and/or 
curvilinear lab frame spacetime coordinates are used, 
and/or the comoving frame undergoes accelerated mo- 
tion, the comoving frame connection coefficients T^gg do 
not vanish, despite the fact that the metric is rj^c, in an 
orthonormal frame. This is because the metric can be 
brought into Minkowski form only locally, with a trans- 
formation that varies from point to point in spacetime. 

Ultimately we want independent variables other than 
the orthonormal comoving frame components x p and p % . 
A global coordinate basis for spacetime position typi- 
cally is needed for numerical simulations, often an Eu- 
lerian one when the particle transport is coupled to mul- 
tidimensional hydrodynamics. As for momentum space, 
Cartesian momentum components normally are replaced 
by spherical-type coordinates, i.e. momentum magnitude 
(or energy, for massless particles) and two angle variables. 

Given a lab frame coordinate basis (denoted by un- 
adorned indices) in which the metric components are 
g^, the coordinate transformation L p p_ to the comoving 
frame is such that 



We note for later use that this implies 
det (L%) = H?)" 1 / 2 , 



(8) 



(9) 



where g — det(g ftiy ). Because the comoving frame 
components of the fluid four-velocity are (u^ — 

. T 



u°, u l ) — (1,0, 0, 0) T , the transformation L p p will in- 



deed be to the comoving frame provided 



L" 8 = lv u 



(10) 



In contrast, L^i are not unique, but are fixed by Eq. ([5]) 
only up to spatial rotations. 

A lab frame coordinate basis also gives us a handle on 
the comoving frame connection coefficients. In the lab 
frame coordinate basis, the connection coefficients are 
torsion- free (symmetric in the lower indices), and given 
in terms of the metric by 



1 



7f 



jtXT 



dg<r V . 9g ap dg up 



dx p dx v 



dx a 



(11) 



The connection coefficients do not transform as tensors, 
but as (see for instance Ref. [2l[ or the Appendix of 
Ref. [H) 

^^i>p = L^^L {,L P p ' V p + pLi p , (12) 

where L p M is the inverse of Note that the comoving 
frame connection coefficients arc not manifestly symmet- 
ric in their lower indices; in comoving frame covariant 
derivatives such as VpU p = dpU^+T^^pU 1 ' , the final index 
of r^j>p should match that of the derivative operators. 

Consider next the use of curvilinear coordinates in mo- 
mentum space (denoted by indices adorned with a tilde), 
for which we use the technology of differential geometry 
in a manner similar to the way it is used in spacetime. 
We introduce a three-metric A^j giving the momentum 
space line element 



d$ 2 = A ? j dp 1 dp 3 



(13) 



where d<& is the 'proper distance' between two points in 
momentum space. The transformation 



dp 1 
dp i 



(14) 



between curvilinear and Cartesian momentum space co- 
ordinates is such that, similar to Eq. (JSJ, 



p l -pi-\ A- 

1 3 ^3 ~ u i3' 



(15) 



where c% is the Kronecker S. Similar to Eqs. (jSJ) and ((9|), 
this implies 



det (P\) = X- 1/2 , 
where A = det (Ajj). The inverse of P l i is 



dp 1 
dp 1 



(16) 



(17) 
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We also define momentum space connection coefficients 
IF-jj,. These can be given in terms of Ajj by an analogue 
of Eq. (fTTj) . Alternatively, an analogue of Eq. (fT2|) gives 



- f)P l ~ 

tji p»_ pj_p^ tt! , pi pk ^ J J 

11 Jfc — r l C ^ k 11 jk * k g j. 

_ pl.pk^ ur 3 
Qpk 



(18) 
(19) 



thanks to the vanishing of the Cartesian connection co- 
efficients IE ~~ k in flat momentum space [HI]. With these 
connection coefficients we define a 'momentum space co- 
variant derivative' operator T>i, giving for instance 



k dp k 3k 



(20) 



for some momentum space three- vector R l . 

We emphasize that the metric Ajj, coordinate transfor- 
mations P l i and P 1 ^, and connection coefficients IPrf in 
momentum space are completely separate from those of 
spacetime. They operate only in (flat) momentum space, 
and we introduce them only to facilitate formal handling 
of generic momentum space curvilinear coordinates. 

As a concrete example of curvilinear momentum space 
coordinates — just an example, meant in no way to spe- 
cialize our formalism to particular momentum space 
coordinates — consider for instance 



p 1 = p cos 

p 2 = p sin $ cos <p, 

p 3 = p sin. i!) sin ip, 



(21) 
(22) 
(23) 



a version of momentum space spherical coordinates. Or- 
dering the curvilinear variables (p*) = (p, ip) T , Eq. (|17|) 
is computed by inspection as 

cos — p sin $ 

[P l i) = | sin $ cos ip poos'!? cos — p sin "8 sin if 
v sin sin ip p cos i9 sin (p p sin i9 cos (p 

(24) 

The transformation P % i of Eq. (fT4")) is then the matrix 
inverse 

I p cos ?9 p sin d cos y p sin d sin y> \ 
(P l i) = — ( — sin?? cos^cost^ cos?9sin(^ 



Finally, 







d$ 2 = dp 2 



— siny/ sin'd cos^/sini^. 



sin 2 -d dp 2 



(25) 



(26) 



is the line element in this example, from which the mo- 
mentum space metric components Ajj can be read. 

With transformations U L ^ and P % i to our desired phase 
space coordinates in hand, we use them to rewrite the 
derivatives in Eq. (JTJ), such that 



Th.ptp^P 1 



' dp 1 



C[f]. (27) 



Thus the dependence of the distribution function is / = 
/ (a; M ,p l ), with lab frame coordinate basis spacetime po- 
sition coordinates x M and comoving frame momentum 
space curvilinear coordinates p % as independent variables. 
The partial derivatives with respect to x 11 are with p l held 
fixed, and vice-versa. 



B. Divergences 

By a 'conservative formulation' we mean one expressed 
in terms of divergences with respect to our chosen phase 
space coordinates. In obtaining a conservative Boltz- 
mann equation we make use of covariant derivatives in 
both spacetime and momentum space. Accordingly, we 
begin by writing Eq. (f2"T|) as 

P^L%W,J - 1% // //'/"/P,/ = C[f], (28) 

using the fact that / is a scalar. Then the product rule 
satisfied by covariant derivatives gives 

-/ V„ +/2* (PT^//) = C[f]. (29) 

The first two terms of this equation are very close to the 
conservative form we seek, and our task is to show how 
the last two 'extra' terms on the left-hand side (almost) 
cancel. 

Consider first the 'extra' spatial derivative term (the 
third on the left-hand side of Eq. ^\)). It is 

I v, (iftpp) = -w ( + r w ) , (so) 



dx^ 



as the partial derivative is taken with p % held fixed. For 
later use we keep in mind that (e.g. Ref. (2lj ) 



vt 1 



1 d. 



-9 



dx v 



(31) 



but we will not need to use this in Eq. (|30j) . 

Processing of the 'extra' momentum derivative term 
fVt {P^iT l p[L p p' 1 )— the fourth on the left-hand side of 
Eq. (|29|) — is interesting but a bit detailed, and we leave 
the specifics to Appendix [A] There are three factors af- 
fected by the derivative. Only in one of them, 2?jP 4 j, do 
the momentum space connection coefficients come into 
play, as this is the only free index (and also the only 
tilde index). It turns out that 



2V\ = o, 



(32) 



which we show in Appendix [Al using Eq. (|19[) . The other 
two factors affected by the derivative are the momenta 
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p v and in 

fp\v- t (r%- lP v) = fP\-^ % (^ 



9p l <9p l 



(33) 
(34) 

(35) 

in which we have used Eq. (|14|) . In the momentum deriva- 
tives, p° must be considered a function of p l through 
Eq. ([3]) . The result for the first term of Eq. ([35]) is worked 
out in Appendix E] (see Eq. (|A28|) ): 



= -/r^ A //p\(-p 8 )A 

The result for the second term of Eq. 
out in Appendix E] (see Eq. (jA34|> ): 



i 



(36) 



3 is also worked 



(37) 



dx^ 



These results are obtained with the help of the geodesic 
Eq. ([6]) , along with a lowered-index version; the fact that 
the mass shell relation of Eq. Q implies = d (p^p* 1 ) = 
2pp_dp tl , where here "d" represents any derivative oper- 
ator; matrix identities involving derivatives, traces, log- 
arithms, and determinants; and Eqs. (O, (O, and ([3T]) . 

We are ready to exhibit our conservative formula- 
tion of the Boltzmann equation, using Eqs. (|30p - ([3"T)) in 
Eq. ([29]). Note that Eq. ([37]) cancels with Eq. ([30]). and 
that Eq. (136)) can be combined with the first momentum 
derivative term in Eq. (|29p . The result is the conservative 
form we seek: 



S N + M N = C[f], 

where 

is the spacetime divergence, and 

1 



M N = (-p 8 ) Vi 



(-po) 



(38) 
(39) 

(40) 



is the momentum space divergence. In terms of partial 
derivatives, the divergences can be expressed 



1 



8 



~N 



M N = 



y/—g dx^ 

(-Po) d 
V\ dp 1 



(-P6, 



(41) 



(42) 



using Eq. pip and its momentum space analogue. In 
comparing with Eq. (132) of Ref. 27], note that (— pg) = 



E(p) and that VA = [det (P\-)] 1 = det (P\) 
t let (dp/du); see Eq. (JTSJ. 

The conservative nature of Eqs. ([55]) . (|4"T]) - (l4*2"]) is seen 
by considering the integral over the phase space volume. 
The invariant momentum space volume element in our 
momentum space curvilinear coordinates is 



dp = 



VA dp 1 dp 2 dp 3 



(-Po) (2^) 3 



(43) 



The factor of vA = det (P\) is the Jacobian of the trans- 
formation from Cartesian momentum coordinates, and 
(— Pq) 1 arises from restricting an invariant 4-momentum 
volume element to the p° > branch of the mass shell 
(e.g. Ref. (3(|). We use units in which % = c = 1; relative 
to works in which instead h — c — 1, this leads to the 
factor of (2tt) 3 in the denominator. When Eq. ([38]) is mul- 
tiplied by dp and integrated over momentum space, the 
momentum space divergence term of Eq. (|42p manifestly 
yields a vanishing surface integral, leaving the number 
balance equation 



1 



d 



-g dx^ 



(V=gN») = J c[f]dp, 



where 



fp^dp 



(44) 



(45) 



is the number flux vector (e.g. Ref. [36j), expressed 
here in terms of the lab frame coordinate basis (note 
pH _ £^ jQ p' 1 ). Turning to position space, recall that 
three- volume elements in spacetime are one-forms whose 
index specifies an orientation. The spatial three- volume 
reckoned by an Eulerian observer (one associated with 
the lab frame coordinate basis) is 



(dx) = y/— g£fj,i23 dx 1 dx 2 dx 3 , 



(46) 



^ vpa is the alternating symbol with £0123 = +1. 
the integrand of Eq. ([45]) , we see that 

(47) 



where e 

Considering the integrand of Eq. (|45|) , we see that 
dN = fp» dp (dx 

VXdp 1 dp 2 dp 3 



(-P )(27t)3 



(y/—g dx 1 dx 2 dx 3 



(48) 



is the number of particles measured by an Eulerian ob- 
server in a differential phase space volume surrounding 
(x l ,p l ) at time t = x . The spatial number density A'" 
of particles measured by an Eulerian observer is 



N°(x*) = / f{x»J) (L%p») 



VX dp 1 dp 2 dp 3 
(-Po)(2^) 3 

and multiplying this spatial density by y/—g dx 1 dx 2 dx 3 
and integrating over all space gives that total number N 
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of particles in the universe represented by the distribution 
function ffx^jp 1 ): 



N(t) 



(IN 



N°(x>*)y/=gdx 1 dx 2 dx 3 



(50) 



Assuming a vanishing particle density at spatial infinity, 
multiplying Eq. (l4"4"]l by ^J—gdx 1 dx 2 dx 3 and integrat- 
ing over all space manifestly yields a vanishing surface 
integral, leaving 



dH(t) 
dt 



= / C[f] 



\f\ dp 1 dp 2 dp 3 



-g dx 1 dx 2 dx 3 



(51) 

for the time rate of change of the total number of particles 
in the universe represented by the distribution function 
fyx^^p 1 ). It is the ready conversion of volume integrals 
to surface integrals via the divergence theorem — allowing 
the time rate of change of a 'conserved' variable in a vol- 
ume to be related to a flux through the surface of the 
volume — that is the hallmark of a conservative formula- 
tion. 



III. 3+1 SPECIALIZATION 

After briefly reviewing the 3+1 formulation of general 
relativity, we introduce 'Eulerian decompositions' of the 
tetrad of comoving frame orthonormal basis vectors; a 
decomposition of the momentum space coordinate trans- 
formation; express the conservative Boltzmann equation 
in terms of the resulting variables; and examine the re- 
lationship between lepton number and four-momentum 
exchange. 



A. Geometry description 

Numerical relativity often is built upon the 3+1 for- 
mulation of general relativity. Pedagogical introductions 
include Refs. HI OH, 53]; see also, very briefly, Sec. IIIA 
ofRef. HI. 

In this approach — which provides our lab frame coor- 
dinate basis — one considers a foliation of spacetime into 
spacelike slices, i.e. three-dimensional hypersurfaces E t 
labeled by coordinate time. The metric components g^ v 
can be read off the line element 

ds 2 = g^ dx^dx" (52) 
= -a 2 dt 2 + 7ij (dx 1 + (3 Z dt) (dx j + /3 j dt) . (53) 

The lapse function a determines the orthogonal proper 
time a dt separating slices E 4 and Ht+dt- The shift vector 
f3 l is tangent to the spacelike slice and characterizes the 
velocity of the spatial coordinates as seen by an observer 
at rest in the slice. The three-metric 7y, also tangent 
to the spacelike slice, encodes its internal geometry. The 
equation 



expresses the determinant g of the four-metric in terms 
of the lapse function and the determinant 7 of the three- 
metric. Another symmetric tensor tangent to the space- 
like slice, the extrinsic curvature K^, describes the warp 
of the spacelike slices as embedded in spacetime. 

In the 3+1 approach, solution of the Einstein equa- 
tions is transformed into a Cauchy problem: specify ini- 
tial data (satisfying certain constraints from the Einstein 
equations) on an initial spacelike slice; and with coordi- 
nate freedom fixed and spatial boundary conditions spec- 
ified, evolve the geometry of the spacelike slices forward 
in time. The dynamical variables are the independent 
components of 7^ and K^, which satisfy first-order- in- 
time evolution equations. The coordinate freedom of 
general relativity corresponds to freedom to specify the 
lapse a and shift p l , which determine the slicing of space- 
time and the motion of spatial coordinates respectively. 
In this paper we regard 7^ and Kij as given, for in- 
stance as having been obtained by numerical solution 
(often of even further transformed systems, as for in- 
stance in Baumgarte-Shapiro-Shibata-Nakamura and re- 
lated approaches); see e.g. Refs. [4CH42I ] . 

'Eulerian decomposition' of tensors into 'Eulerian pro- 
jections' orthogonal and tangent to the spacelike slice will 
prove useful. (This served us well in Ref. [35| on a mo- 
ments approach to radiation transport.) The unit normal 
n M to a spacelike slice — which is also the four- velocity of 
an 'Eulerian observer' — has lab frame coordinate basis 
components 



K) 



(l/a,-/3'Va) T , 
(-a, 0,0,0). 



The orthogonal projector is 



-g = a yfy 



(54) 



(55) 
(56) 



(57) 



From Eqs. (|53|) and (|56j) it follows that the spatial part 
of 7^,^ equals the three- metric 7^ , motivating use of the 
same base symbol. Contraction of an arbitrary vector 
with n M yields the portion orthogonal to a spacelike slice, 
and contraction with 7/i „ = g^ + n^n^ yields the portion 
tangent to the spacelike slice. Indeed a trivial calculation 
confirms that "f^n u = 0. 



B. Tetrad decomposition 

In Sec. [IT] we mostly thought of L^p, as a transforma- 
tion matrix between the lab frame coordinate basis and 
the orthonormal comoving basis, but we can focus in- 
stead on the fact that it represents the tetrad of comov- 
ing orthonormal basis vectors. In this perspective, the 
upper index p, in L 11 ^ labels the (lab frame coordinate ba- 
sis) component of the comoving orthonormal basis vector, 
while the lower index p, is not a component but simply a 
label identifying the particular vector in the tetrad. 

We already saw something of this in Eq. (fTU)) . and we 
can extend the thinking to the comoving spatial unit vec- 
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tors. We saw in Eq. (fT0|) that the components of the La- 
grangian observer's four-velocity in the comoving frame 
are 

(^) =(1,0,0, 0) T , (58) 

so that 

L" 6 = L" d u 6 = L^u* = (59) 

That is, L^g is the Lagrangian observer's timelike basis 
vector, expressed in components measured by an Eulerian 
observer. Similarly, the Lagrangian observer can choose 



spatial basis vectors 

(w*) = (0,1,0, Of, (60) 

(/) = (0,0,1, Of , (61) 

= (0,0,0, if , (62) 
whose lab frame coordinate basis components are 

L^j = L^iw' 1 = = w^, (63) 

= lV = l^v* = ( 64 ) 

L^3 = L"§2 S = = z». (65) 



For the time being we will use the individual expressions 
m m , z M to emphasize that we are working with a 

tetrad of vectors. Later, for the sake of more compact ex- 
pressions using index notation, we will return to notation 
like L^p,, and ask that the reader keep expressions like 
Eqs. ([5^1 . (|6"3"l) - ([6"5)) in mind, and regard hatted indices 
like jl only as labels identifying a particular vector in the 
tetrad in such instances. 

In a moments approach to radiation transport [35| . 
the only comoving basis vector we had to deal with was 
L m q = u M , and we found it useful to consider its Eulerian 
decomposition into parts orthogonal and tangent to the 
spacelike slice: 

u" =A(n" + u") . (66) 

The orthogonality requirement on v^, 

= 0, (67) 

implies (see Eq. (|5d| ) that is spacelike and has com- 
ponents 

K) = (0Xf (68) 

in the lab frame coordinate basis. The interpretation of 
as the three- velocity of a Lagrangian observer as mea- 
sured by an Eulerian observer is confirmed by squaring 
Eq. (|6"6")l to find the expected Lorentz factor 

A = (1 - «%r 1/2 = (1 - v* Vi y 1/2 (69) 

for a boost between frames with relative three-velocity 
(Recall that u^u^ — n^n^ = —1.) 



Here we also find it useful to make Eulerian decompo- 
sitions of the spatial comoving basis vectors u^, y M , z M . 
We write 

= An^ + , (70) 
y" = Bn" + b", (71) 
z^ = CV l + c^, (72) 

specifying that 

= n^aF- = = n^, (73) 

or, in the lab frame coordinate basis, 

K) = (OXf, (&**) = (O^f, (c») = (0,c l ) T (74) 

in order that a M , b^, d 1 be tangent to the spacelike slice. 

The metric transformation in Eq. (|SJ) translates into or- 
thonormality conditions on it M , w^, y^jZ^. As a physical 
four- velocity, the timelike basis vector u M already satisfies 

= L^L\ 9fiV = m = -1. (75) 

The off-diagonal time-space elements of Eq. © , 

L^ b L v t g lu> = ri 6t = 0, (76) 

give orthogonality relations 

= = = u^, (77) 

which turn out to imply 

A = « M a" = v,a\ (78) 
B = = Vib\ (79) 
C = v„cf = v iC \ (80) 

where we have used Eq. (j6"B"|) . The diagonal spatial ele- 
ments of Eq. ([8]) , 

L^iVi = rju = +1 (no summation), (81) 

impose the normalization conditions 

+ 1 = Wfi w» = y^ = z^, (82) 

which imply 

lij a l a 3 '' = cua 1 = a^a? = 1 + A 2 , (83) 
7y = bib* = bfJT = 1 + B 2 , (84) 
jij cV = ad - c„c" = l + C 2 (85) 

for the normalization of 0,^,6^,0^ (the portions of 
, y M , z M tangent to the spacelike slice). Finally, the 
off-diagonal spatial elements of Eq. ((SJ, 

L^ 3 g^= ms = (i^j), (86) 

give orthogonality relations 

= = = y^, (87) 
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which imply 

jij a l V = a t b l = = AB, (88) 

7ij dV = a lC l = a M c M = AC, (89) 

jij 6V' = & lC * = &„c" = SC. (90) 

Thus orthogonality of w^^^^z^ requires that their 
parts a^,b^,c^ tangent to the spacelike slice are not 
orthogonal — for v % ^ 0, that is. For v % — > 0, the orthonor- 
mal comoving tetrad u 11 , , y 11 , z^ becomes an orthonor- 
mal lab frame tetrad: the timelike basis vector — > 
(Eqs. (|66|) . ([69])). becoming orthogonal to the spacelike 
slice; and A, B, C ->• (Eqs. ([78]) - ([SD ]) ). so that the or- 
thonormal triad w^,y^,z^ a^,b^,c^ (Eqs. (|70>([72|) ). 
with the tangents a' 1 , & M , c^ 1 to the spacelike slice becom- 
ing orthonormal (Eqs. ([55 ]) -(|55 ]) , dHHJ) - dHOJ) ) ■ (Of course, 
this v l —> lab frame orthonormal basis generally does 
not coincide with the non-orthonormal lab frame coordi- 
nate basis, as is evident from the fact that Eq. (|53[) is not 
the Minkowski metric for generic g^ v .) 

Specification of the comoving frame tetrad 
vf* , , y^ , is not yet complete. The timelike 
comoving basis vector — the Lagrangian observer four- 
velocity — is determined by the evolution of the fluid 
and geometry. As for the spatial comoving basis vectors 
it^, z M , we have given their Eulerian decompositions 
in Eqs. ([70]) - (172"]) . We can regard A, B, C as determined 
by Eqs. ([75 ]) -([SO ]) , In light of Eq. j7J|, this leaves nine 
unknowns — a 1 , b l , c l — but only six equations so far for 
them (Eqs. ([85 ]) -([85 ]) , ©- flM])). The remaining three 
degrees of freedom correspond to the choice of spatial 
orientation of the comoving orthonormal frame. 

One possibility for pinning down the spatial orientation 
would be to insist that the Lagrangian observer Fermi- 
Walker transport her orthonormal triad of spatial ba- 
sis vectors, disallowing spatial rotations in the comoving 
frame [2l[. This simplifies life for the Lagrangian ob- 
server by removing Coriolis force contributions to the co- 
moving frame connection coefficients T^cp. It would also 
provide a means of eliminating the appearance of the time 
derivatives of the tetrads we shall encounter in the mo- 
mentum space divergence in Eq. (|4"2"]l . However, it would 
add to the list of evolution equations to be solved, and 
complicate the connection between the particle number 
and four-momentum conservative reformulations of the 
Boltzmann equation (see Sec. IIII E]) . 

We choose instead to specify three algebraic conditions 
that simplify the diagonalization of jij implicit in Eq. ([5]) . 
It could be required, for instance, that certain lab frame 
coordinate basis components of b l and c* vanish: 

6 X =0, (91) 
c 1 = c 2 = 0. (92) 

The remaining components of a % ,b l , c l are then deter- 
mined by (a) a single equation to be solved for the single 
non- vanishing component c 3 of c l ; (b) two equations to 
be solved for the components b 2 and b 3 of b l ; and (c) three 
equations to be solved for all three components of a 1 . The 



resulting equations are not particularly illuminating, and 
we relegate them to Appendix [B] We do note that in the 
special case of diagonal and v l = 0, these equations re- 
duce as expected to a 1 — (7n) _1//2 , b 2 = (722) 1//2 , c 3 = 

j/2 

(733) i with all other components vanishing. 

For use in more compact index notation, it is conve- 
nient to collect the Eulerian decompositions in Eqs. (|66p . 

as 

L" A = n%i + .P A , (93) 

where 

(C (i ) = (A,A,B,C), (94) 
= (Aif^V,*), (95) 

keeping in mind that t° p, = according to Eqs. ([68]) and 
([74]) . The choices in Eqs. (HJ) and |H2]) make £*| (viewed 
as a matrix) lower triangular. (We caution that the use 
of the base symbol I in Ref. [35[ differs from its use here.) 
Having discussed the interpretation of 

(L^) = (u»,w»,y»,z> 1 ) (96) 

as a tetrad of comoving frame orthonormal basis vectors, 
it is easy to see that 

(L^p) = {-Up, Wp,yp, Zp) T (97) 

is the inverse, for the orthonormality conditions render 
Eqs. ([M]) and (H7J consistent with 

L%L^ = 5^. (98) 

The Eulerian decomposition of the inverse is 

L\ = £»np + £^, (99) 

with 

{&) = {-K,A,B,C) T , (100) 
(£%) = (-Avp,ap,bp lC pf . (101) 
In addition to Eq. ([9"B"]) . we have also 

L^lK = 5» v . (102) 

Inserting Eqs. ([Ml) and in Eq. ([TU2"]) and taking ap- 
propriate contractions with n M and ~ip V , we find 

= -1, (103) 

C^p = 0, (104) 

i 11 ^ = 0, (105) 

*V A * = 7"* (106) 

for various contractions of the Eulerian projections of 
L^p and L^p. These expressions of tetrad orthonormal- 
ity in terms of the Eulerian projections will prove useful 
in Sec. iHIEl 
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C. Momentum decomposition 

In a spirit similar to the decomposition of TJ 1 ^ — and for 
reasons that will be become more apparent in Sees. 1111 Dl 
and MI El — we will hnd it useful to decompose the mo- 
mentum space transformation P 1 \ and its inverse P l i (see 
Eqs. (|14p and (fTT)) ) into pieces parallel and orthogonal to 
the particle three- momentum p 1 . In particular we write 



P l , = 



IP-,. 



Here 



P 

is the projection of P % i parallel to p l , and 



P 



<p l Pi 



(107) 



(108) 



(109) 



is the three-momentum magnitude. The projection or- 
thogonal to p % is 



U % P jk^ij 



(110) 



^ given in Eqs. (l93l) and (|9Q|) . These appear in the 
spacetime divergence in Eq. (|41[) and in the comoving 
frame connection coefficients in the momentum space di- 
vergence in Eq. ((32). The expression for the comoving 
frame connection coefficients in Eq. (fl~2j) can be recast as 



>'i> 



(118) 



This is the expression for the 'Ricci rotation coefficients' 
given by Lindquist [l(| ■ Riffert [IH emphasized Eq. (fl"2"|) , 
claiming greater simplicity, but their equivalence is ap- 
parent. For our purposes, Eq. (I118[) is more conve- 
nient, as it allows us to obtain a 3+1 formulation of 
the Boltzmann equation while almost completely avoid- 
ing explicit encounters with connection coefficients (see 
Appendix [C| . 

Of the two terms on the left-hand side of Eq. (1381) . 
the spacetime divergence given by Eq. flU} is by far the 
simpler. Using Eqs. (j5~4"l) . (|9"3")l , and ([55]). this spacetime 
divergence is 



S 



N 



("Pel 



o(d n ) , d(F N y 



dt 



dx l 



(119) 



in which 



where 



P°Pi 



(111) 



is a momentum space orthogonal projector. Similar to 
Eq. (|107|) . the decomposition of the inverse is 



(FnY 



V7 



(-P6) ~" 
(-Po) 



CtfPf, (120) 
(atfi-FCdpPf (121) 



P\ = 



p'Qi 
p 



(112) 



As with and its projections in Eqs. (|102p - (|106[) . the 
fact that 



(113) 



implies the following relations among the projections: 

QiQ 1 = 1, (H4) 



U l iQ l = 0, 
Q%U\ = 0, 
U %U j = k j, 



(115) 
(116) 
(117) 



which follow from inserting Eqs. (|107l) and (|112j) into 
Eq. (|113|) and taking appropriate contractions with p l 



and fc l j. 



D. Boltzmann equation 

We now give, in 3+1 form, the divergences in the 
number-conservative formulation of the Boltzmann equa- 
tion in Eq. (|3"8"|) . Key to our derivation is use of the Eu- 
lerian decompositions of the tetrad p, and its inverse 



are respectively the conserved number density and num- 
ber flux. 

Using the results for the comoving frame connection 
coefficients given in Appendix [C] the momentum space 
divergence in Eq. (|42p can be expressed 



M 



N 



1 (-Po) 9 
"V7 VX dp 1 



41- 



'(-Pd) 



(Rn) 



{R 



N 



(o N y}}, 



(On, 
(122) 



where 
(RnY 



P jrf 



(-Po) 



C p P ( 

a dxi 



l k ~* K jk 



CoCp, da e k u£p,d(3 k tot^d-fk, 



dxi 
(123) 



describes momentum shifts (that is, redshift and angu- 
lar aberration in momentum space spherical coordinates) 
due to gravity as embodied in the spacetime geometry, 
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and 



(On) 



p _ 



V7 

(-pa) 
c p 



v v v p f 



C 



dt 
dt 



+ 



(at 3 a 



^) ft* 



dC 

de kl 



(124) 



are 'observer corrections' due to the acceleration of a La- 
grangian observer (and partially entangled with the ge- 
ometry as well), arising from the tetrad connecting 
the orthonormal comoving frame to the lab frame coor- 
dinate basis. It is through the decomposition of P l % in 
Eq. (|T07l) — along with Eq. (|C24|) — -that the hrst term of 

Eq. (I122p contains r"j>£, expressed here in terms of (Rn)° 

and (O n ) v . This will prove beneficial in relating the 
number-conservative formulation of the Boltzmann equa- 
tion to the four-momentum-conservative formulation. 



E. Lepton number and four-momentum 
reconciliation 



In addition to the number-conservative reformulation 
of the Boltzmann equation given by Eqs. (|38l) - (|4"2"|) . there 
is also a four-momentum-conservative reformulation [27| . 
It has the same basic structure, with spacetime and mo- 
mentum space divergences, but with an additional index: 



(s T y + (m t ) p = p pc[f] 



(125) 



Each of the components of this equation is redundant 
with Eq. (l38l) for number exchange with the fluid; but 
the right-hand side yields source terms for the fluid 
equations, and the structure of the left-hand side illu- 
minates issues faced in accounting for energy and mo- 
mentum conservation in the discretized case. We ex- 
amine the forms taken by (Sr) p and (MtY in our 
chosen phase space coordinates, and consider how the 
four-momentum-conservative reformuation in Eq. (|125[) 
is related to the number-conservative reformulation in 
Eq. O. 

We begin at a high level, working with covariant 
derivatives in spacetime. Comparing the right-hand sides 
of Eqs. (|125l) and Eq. (f38|). we see that the left-hand side 
of Eq. (I125[) must follow from multiplying the left-hand 
side of Eq. (|3"51) by p p = L p p p p . Using Sn as given in 
Eq. flSH]), we have 



where 



and 



p p S N = (S T ) P + (Es) f 



(5t) p = V m (L^L^p^f), 



(Es 



-fL^ptpP (V^Z/p) 



(126) 



(127) 



(128) 



is the 'extra' term that results from pulling p p inside the 
divergence. Turning to the momentum space divergence 
and using Mjy as given in Eq. (|4*2"j) . we have 



p p M N = (M T ) P + (EmY 



(129) 



where 



(M T y = 



(-Po) d 



(-Po) 



PhT^^p^LPppPf 



and 



(E 



M 



dp p 



dpP 
dp 1 



— fV l ~~ ■n v rJ i T P ~ — 

- fl fyP p L, p -g-j, 



(130) 

(131) 
(132) 



in which we have used Eq. (|14[) in the second line. Using 
Eqs. (lAip . (|A13[) . and (|6|), this 'extra' momentum space 
derivative term becomes 



This can be expressed as 

(E M Y = fL^ (VpL p )p ,y p^ 



(133) 



(134) 



by virtue of Eqs. (fill 
Eq. (fT28l) . we see that 



and (|102[) . Comparing with 



(EsY + (E 



M 



0. 



(135) 



Thus Eq. (|125p is verified with the divergences (StY an d 
(M T y given by Eqs. (fHT)) and (fT3D|) respectively. 

The relationship between Eq. (|38[) for particle number 
exchange with the fluid and Eq. (j 125[) for four-momentum 
exchange also can be examined in terms of the more 
detailed expressions for Sn and Mn we have given in 
Eqs. (1119p - (1124[) . In doing so it is convenient to split the 
four-momentum equations into an energy equation and a 
three-momentum equation. The equation for energy ex- 
change with the fluid is obtained by multiplying Eq. (|3"5)) 
by the lab frame neutrino energy, i.e. the projection of 
p p orthogonal to the spacelike slice: 



n p p F 



i p L P pp p = £p P p , 



(136) 



where we have used the Eulerian decomposition of L p p 
given by Eq. (|93|) . An equation for momentum exchange 
with the fluid is obtained by multiplying Eq. (|38| by the 
lab frame neutrino three-momentum, i.e. the projection 
of p p tangent to the spacelike slice: 



l3 P P P = l 1P L%P p =t 3P P p , 



(137) 



where once again we have used the Eulerian decompo- 
sition of L p p given by Eq. (|93j) . Let us represent the 
right-hand sides of both Eq. ([He"]) and (fT3?|) by 

£„,/•. / ,„/-" ■ /.pif. (138) 
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Thus using Zp p p to represent both the lab frame energy 
CpP p and the lab frame momentum £jpp p , we derive in 
parallel the energy- and momentum-conservative refor- 
mulations of the Boltzmann equation. 

We begin with the spacetime divergence. First, using 
Eq. (|TT9 j) -(pI ]l . we have 



Zp P p S_ 



N 



(-po) 



T.Z 



T.Z) 



dt 



dx l 



+ Es,z 



Here 



D T ,z = Z p pPD 



N 



(139) 



(140) 



is the 'conserved' density; for Z p p p — > C p p p we have the 
'conserved' energy density 



D T ,n = £pP p D N , 



(141) 



Next we turn to the momentum space divergence. The 
result of multiplying Eq. (|122[) with Z^p p can be cast in 
the form 



Zpp p M. 



N 



Here 



i (-go) 



+ 



(-Po) 



P 

(Rt,z) 



(Rr,z)° + (Ot,z)° 



(0 T ,zf\ } 



Em,z- 



(R T ,zf = Z- p pP {R N f . 
(0 T , z f = Zp P p (0 N f. 



(149) 

(150) 
(151) 



while for Zpp p — ¥ ljpP p we have the 'conserved' momen- 
tum density 



., £jpP p D N . 
Next in Eq. (fl3"9"j) is the flux 



{F TiZ y = Z pI f(F N y, 
which for Zpp p — > CpP p is the energy flux 

(F T , n )* = Cpp? (F N )\ 
and for Zp p p — > ijp p p is the momentum flux 



(142) 



(143) 



(144) 



(145) 



Finally in Eq. (|139[) we have 



E. 



s,z 



(-Po) 



p p p p f 



dZ a 



(146) 

the 'extra' term that comes from pulling Zp p p inside the 
time and space derivatives. For Z p p p — > CpP p , 



Es,n 



(-Pa) 



r d S± 
dt 



(at 



(147) 

is the 'extra' term in the energy equation, while for 

ZppP^ljppP 



(E Srt ) 



V7 

(-Po) 
d^p 

p dt 



p p p^f 



(148) 



is the 'extra' term in the momentum equation. (We cau- 
tion that E s , n ^ -n p (E s ) p and (Es^j ^ -j jp (E s ) p , 
where {Eg) p is given by Eq. (|128[) . as might plausibly but 
incorrectly be suggested by the notation.) 



(R T ,nY =C pP p (R N 



p- - r ,,/' i v) . i'' 
p- — a*, (n..\v 



(Ot,v) — £pP p (On) 
for Zpp p — » £ ( sp p (energy equation), and 

(RTnfj =ly p p p (R N f, 



W (Oat) 



(152) 
(153) 

(154) 
(155) 



for ZpP p — > ^jpP p (momentum equation), with (Rn) p 
and (O n ) p given by Eqs. (fT2"3"|) and (fT^)> respectively. 
The 'extra' term E M z in Eq. (|149[) is 



Q ? (-Po) 



(iijv) + (Oat) 



(R N y + (o N y]}z p ^. (156) 



The last two factors can be rewritten as 

dp p _ P j 7 dp p 
p d P > ~ dpi 



P>Qi 



Zn 



Pj 



(-Pa) 



(157) 
(158) 



in which we have used Eqs. JT7|}, ([IT2]) , and (TXTJ. Plug- 
ging back into Eq. (fTM)) and using Eqs. (fIT4| - (fTT7l) yields 



E 



M.Z 



(-Po) P 



3 r 



P P 



(R N )° + (O n )° 



(R N y + (On)"] } 



Zf, 



Pj 



°(-P6) 



Zi 



(159) 



This reduces further to 
Em,z = 



y P'Zj 



(R n ) p + (O n ) p (160) 



= -z P 



(R N y + (o N y 



(161) 
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The second line follows from the fact that 



o = P pr%pV\ (162) 



by virtue of Eqs. (|A13|) and ([5]), through which 



0=Pp (R 



IN 



(OnY 



(163) 



is also implied (see Appendix [Q and Sec. IIIIDj) . 

Considering the first term in Eq. (|161l) , we will see that 
the contribution of the redshift/aberration terms (Rn) p 
to the 'extra' term Em,z from the momentum space 
divergence becomes the geometric/gravitational source 
term arising from the covariant derivative in Eq. (|127l) . 
For Zp — > Cp (energy equation), we have 



Gr,n — —£p {Rn) p 



{-Po) \ a 0x3 



(164) 
-l k pK jk ){\Qh) 



in which only the second line of Eq. (|123p contributes, 
thanks to Eqs. (|T0"3"|) and (fl"04|) . On the other hand, for 
Ijp (momentum equation) we have 



(Gt, 7 ) 7 - - ~ijp (RnY 



(-Po) 



df3 k 



LpCp da 

a dxi 

l k o Pp. dj ki 
2 dxi 



(166) 



(167) 



in which only the third line of Eq. (II 23|) contributes, 
thanks to Eqs. (fTUS) and (fT06]> . 

Turning to the second term of Eq. (|16ip , we show that 
the contribution of the observer corrections (Om) p to the 
'extra' term Em.z from the momentum space divergence 
cancels with the 'extra' term Es,z from the spacetime 
divergence. Indeed, for Zp — > Cp (energy equation), 
only the second line of Eq. (|124j) contributes thanks to 
Eqs. fT03]) and (|T0l]> : 



— Zp [OnY — —Es, n , 



(168) 



where we have compared the result with Eq. (|147j) . On 
the other hand, for Zp —> £jp (momentum equation), 
only the third line of Eq. (|124j) contributes thanks to 
Eqs. ([T05]) and (|T06| : 



t iP (OnY 



(169) 



where we have compared the result with Eq. (|148p . 

We are now in a position to combine the above results 
into the energy- and momentum-conservative reformula- 
tions of the left-hand side of the Boltzmann equation. In 
the energy equation, 



Cp P p (S N + M N ) 



-n pP P(S N + M N ) (170) 
-n p [{S T Y + (M T Y\ , (171) 



where 

- n p (S T ) 



, _ (-Pa) 



d(Ih 



d(F 7 



dt 



dx l 



and 



Gr,n 
(172) 



(R T ,n)° + (0 T ,nY 



-n p (M T Y = ^- ( -^- 
VA 

dp 1 [ p 

+V\U\ [{R T , n y + (Or,n)*] } ■ (173) 

We note the consistency of Eq. (I172p with Eq. (A14) of 
Ref. HI, with the substitution -> pL v p il p £ '} '. In 
the momentum equation, we have 

tjfi P p (Sn + M N ) = 1]P p p (S N + M N ) (174) 
= 1 JP [(S T Y + (M T Y], (175) 



where 



Hp(St) p = 



-Pa j 



*V7 



dt 



dx l 



and 



- (Gt, 7 ) 3 

(176) 



hp (M T y 



1 (-Po) 
op 



f(-Pd) 



(RtJ 3 + (Qr, 7 )° . 



+V\U\ [(Rr^Tj + (OtJ] } • (177) 

We note the consistency of Eq. ([776]) with Eq. (A31) of 
Ref. [H, with the substitution T» v -> L p op^p f . 



IV. CONCLUSION 

We conclude by assembling expressions obtained in 
Sec. IIIII into conservative 3+1 general relativistic refor- 
mulations of the Boltzmann equation, and discussing lep- 
ton number and four-momentum exchange between neu- 
trinos and the fluid. Tables Hl lIVI present overviews of 
many of the variables that have been assembled into the 
major entities appearing in these equations. 

The number-conservative 3+1 general relativistic 
Boltzmann equation is 



o(d n ) , d(F N y 



dt 



dx" 1 
1 d 



(-Po) 



(R N )° + (0 N f 



\f\dp l 

+Vxu'\ [(R N y + (o N y] } = c N . (ire) 

This is a conservation law (or balance equation) for par- 
ticle number. The particle density Djq and spatial flux 
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(F N f are given by Eqs. fjl20[) and (|121|) respectively. The 

contributions (Rn) P and (Oat) p to the momentum space 
flux — which represent momentum shifts (redshift and an- 
gular aberration) due to geometry /gravitation and accel- 
eration of a Lagrangian observer respectively — are given 
by Eqs. (fT23]> and (fT24]) . The right-hand side is 

C N = p^C[f], (179) 

where a and 7 are respectively the lapse function and the 
determinant of the three-metric characterizing the inter- 
nal geometry of a spacelike slice (see Sec. IIII A[) . (— pg) 
is the particle energy in the comoving frame, and C[f] is 
the invariant collision integral appearing for instance in 
Eq. (0). 

The phase space coordinates used in Eq. f| 1 78[) al- 
low particle/fluid interactions to be evaluated in the co- 
moving frame in the context of Eulerian grid-based ap- 
proaches to multidimensional spatial dependence. The 
global 'lab frame' spacetime coordinates t and x % are 
those associated with the 3+1 formulation of general rel- 
ativity, in which the line element and metric components 
g^u are given by Eq. (|53l) . In momentum space we ap- 
ply the transformation P l i — dp 1 /dp 1 from comoving 
frame orthonormal Cartesian momentum components p l 
to comoving frame curvilinear momentum components 
p l . This transformation is decomposed into a portion Q % 
parallel to p l , and a portion U % % orthogonal to p l ; see 
Eq. (|107[) . For the case of massless particles, we note 
that the particle energy (— P5) cancels with the three- 
momentum magnitude p = y 1 pip 1 in Eq. (|178[) . In our 
derivation we introduce a momentum space metric with 
determinant A (see the 'line element' in Eq. (fTB"]) . and 
also Eq. (1161) ). A particular example of curvilinear mo- 
mentum space coordinates is the version of momentum 
space spherical coordinates discussed in the next-to-last 
paragraph of Sec. Ill Al 

Making the Boltzmann equation explicit has required 
specification of L^p,, which can be regarded either as the 



TABLE I. Some spacetime and fluid variables. 


E t 


Spacelike slice labeled by time t 


Sec.|IIIA| 


x» 


Spacetime position 


Sec.|IIlA| 


a 


Lapse function 


Eq. §3} 


P* 


Shift vector 


Eq. J53} 




Three-metric 


Eq. (EH) 


7 


Determinant of 7^ 


Eq. (SHI 


Kii 


Extrinsic curvature 


Eq. ([C2| 


n" 


Unit normal to Et; Eulerian four- velocity 


Eq. (H5j 




Lagrangian (i.e. fluid) four-velocity 


Eq. dm) 




Three- velocity of Lagrangian observers 


Eq. J66]) 


A 


Lorentz factor of Lagrangian observers 


Eq. J69l) 


7m- 


Projector orthogonal to n M 


Eq. (57]) 



transformation from a comoving frame orthonormal ba- 
sis to the lab frame coordinate basis, or as the tetrad of 
comoving frame basis vectors; see Sec. IIII Bl The timelike 
vector L^q = is uniquely specified as the four-velocity 
of a Lagrangian observer riding along with the fluid, but 
additional conditions must be imposed to fix the orien- 
tation of the triad of comoving spatial basis vectors; we 
have chosen algebraic conditions that simplify the trans- 
formation to an orthonormal frame. Also, in order to 
arrive at a useful 3+1 formulation we have found it con- 
venient to perform an 'Eulerian decomposition' of L^^ 
into its 'Eulerian projections' and £ l p, orthogonal and 
tangent to the spacelike slice respectively. This geometric 
approach (in perspective if not notation) enables evalu- 
ation of the momentum space divergence while almost 
completely avoiding explicit encounters with connection 
coefficients. 

Turning to astrophysical simulations in which neu- 
trinos interact with matter, we note that solution of 
Eq. (|178p — for Dn, and ultimately for the distribution 
function / (a:^, p 1 ) , through Eq. (|120[) — enables calcula- 
tion of lepton number exchange with the fluid. Conser- 
vation of total electron lepton number is expressed 

V M (JV£ + N* - JV&) = 0, (180) 

where N£ is the net electron number flux vector of the 
fluid, and N£ and are the number flux vectors of the 
electron neutrinos and antineutrinos respectively. From 



TABLE II. Transformation L M ^ between comoving and lab 
frames, which also can be regarded as the tetrad of comov- 
ing orthonormal basis vectors; and its Eulerian decomposition 
into projections normal and tangent to the spacelike slice E t . 



7"M . 

Lt p. 


Transformation, or comoving tetrad 


Eq. 


|96l) 




L M q; timelike comoving unit vector 


Eq. 


(59]) 




L M j; spacelike comoving unit vector 


Eq. 


(|63|1 


v" 


L^2\ spacelike comoving unit vector 


Eq. 


([64]) 




L M 3; spacelike comoving unit vector 


Eq. 


(ED) 


Li jx 


Eulerian decomposition 


Eq. 


([93]) 


u» 


Eulerian decomposition 


Eq. 


dSSJ) 




Eulerian decomposition 


Eq. 


Cm]) 




Eulerian decomposition 


Eq. 


CD 




Eulerian decomposition 


Eq. 


(O 




Projection of L M p normal to Et 


Eq. 


(El 


/>M . 

1 M 


Projection of L M /i tangent to Et 


Eq. 


(ESI) 


A 


Projection of u M normal to Et 


Eq. 


PI 


A 


Projection of ui M normal to Et 


Eq. 


dZSJ) 


B 


Projection of j/ M normal to Et 


Eq. 


US]) 


C 


Projection of z M normal to Et 


Eq. 


Ho} 


Av l 


Projection of it M tangent to Et 


Eq. 


(E3) 


a* 


Projection of w M tangent to E t 


Appendix [Bl 


V 


Projection of y M tangent to Et 


Appendix |B1 


e 


Projection of z M tangent to E t 


Appendix IBI 
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Eq. (|180[) and the discussion in Sec. Ill Bt we see that 
the fluid electron lepton number source term is the mo- 
mentum space integral of the difference of the v e and v e 
collision integrals: 

V„AT = - J {C Ve [/] - C Pc [/]) dp. (181) 

Because the left-hand side takes the conservative form 
1 d 



(182) 



solving for the neutrino distributions using Eq. (|1T8|) 
which is manifestly conservative in both position 
space and momentum space — renders almost trivial the 
achievement of machine-precision global lepton number 
conservation (the position space integral of Eq. (11801) ) in 
numerical simulations. 

In addition to the number-conservative formulation in 
Eq. (|178l) . we have also the energy- and momentum- 
conservative formulations 



d(D 7 



d(F T 



Of 



J__d_ 

= Gt.u + Ct,ui 

d{D Tn ) j d(F T „y. 



(R T , n )° + (0 T , n )° 



(Or J}} 



(183) 



Ot 



dx l 



(i?T, 7 r>(o T , 7 ) 5 j} 



(G T , 7 ), + (C r , 7 ), 



(184) 



TABLE III. Particle momentum; particle distribution func- 
tion; transformation to comoving curvilinear momentum co- 
ordinates, and its decomposition and projections; curvilinear 
momentum space metric. 





Coordinate basis 


Sec 


|IIA| 


P A 


Comoving orthonormal basis 


Sec 


|IIA| 


p l 


Comoving curvilinear basis 


Sec 


|IIA| 


V 


Three-momentum magnitude 


Eq. 


(fT09} 




Particle distribution function 


Sec 


|IIA| 


PS 


Transformation to curvilinear basis 


Eq. 


m 


P\ 


Decomposition with respect to p 1 


Eq. 




Q l 


Projection of P l t parallel to p 1 


Eq. 


(fi"08} 


U\ 


Projection of P % % orthogonal to p % 


Eq. 


(fTTU} 




Projector orthogonal to p 1 


Eq. 


CD 




Curvilinear momentum space metric 


Eq. 




A 


Determinant of \ij 


Eq. 





The energy and momentum densities (Dt,u) and (Dx )7 )j 
are given by Eqs. (|141[) and (|142[) respectively, and the 
corresponding spatial fluxes (Ft.tiY and (Ft^Y ■ are 
given by Eqs. (|144|) and (|145j) . The contributions {Rt.uY 
and {Rt.jY • to the momentum space flux — which repre- 
sent momentum shifts (redshift and angular aberration) 
due to geometry/gravitation — are given by Eqs. (|152[) 
and (|154|) respectively; and the contributions (Ot,«) p and 
{Qt,iY j ^ ue t° the acceleration of a Lagrangian observer 
are given by Eqs. (|153[) and (|155[) . The source terms 



Cn 



(-Po) 



(C , T, 7 ) j = r~\ P^ c lf] 
\ Po) 



(185) 
(186) 



are due to particle interactions. The energy- and 
momentum-conservative formulations in Eqs. (|183l) and 
(|184p contain extra source terms Gt,u and (C?t, 7 )-, 
given by Eqs. (|165[) and (I167[) , relative to the number- 
conservative formulation in Eq. (|178[) . Physically, these 
source terms represent the exchange of energy and mo- 
mentum with the gravitational field; mathematically, 
they are a consequence of the fact that Eqs. (|183[) and 
(jl 84[) are projections of the divergence of a rank-two ten- 
sor, which produces an additional connection coefficient 
term relative to the divergence of a vector. (In our deriva- 
tion in Sec. lIIlEl we do not compute these connection co- 
efficients directly; instead, Gt,™ and (Gr, 7 )j emerge from 



TABLE IV. Major entities appearing in Eq. JT78J, the 3+1 
general relativistic number-conservative reformulation of the 
Boltzmann equation; and in Eqs. (|183p and (|184[1 . the 3+1 
general relativistic energy- and momentum-conservative re- 
formulations of the Boltzmann equation. 



D N 


Conserved number density 


Eq. 


(fT20l) 


D T ,„ 


Conserved energy density 


Eq. 


(fTiTj) 


{Dt, 1 ) j 


Conserved momentum density 


Eq. 




(FnY 


Number flux 


Eq. 


(fmj) 


(FT,nY 


Energy flux 


Eq. 


cup 


(Ft,,)^ 


Momentum flux 


Eq. 


CHS} 


(RnY 


Gravitational shifts, number 


Eq. 


GUI) 


(Rt.uY 


Gravitational shifts, energy 


Eq. 


CE52J) 


{Rt,iY j 


Gravitational shifts, momentum 


Eq. 


CESH 


(OnY 


Observer corrections, number 


Eq. 


ipl} 


(Ot.uY 


Observer corrections, energy 


Eq. 


(fT53)l 




Observer corrections, momentum 


Eq. 


CESS} 


Gt,ti 


Gravitational energy source 


Eq. 


cess} 


(Gt„Y 


Gravitational momentum source 


Eq. 


(fl67|) 


Cn 


Collision number source 


Eq. 


CEZU) 


Ct,ti 


Collision energy source 


Eq. 


CESS} 




Collision momentum source 


Eq. 


CEM} 
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the second and third lines of (Rn) p in Eq. (|123[) respec- 
tively when the momentum projections p^ and £jn p^ 
orthogonal and tangent to the spacelike slice — which re- 
late Eqs. ([155} and (fI5S|l to Eq. ([175} — are pulled into 
the momentum divergence.) 

Consideration of the particle stress-energy tensor 
shows how the Boltzmann equation is related to the ex- 
change of energy and momentum with the fluid. While 
Eq. (|45[) gives the particle number flux vector in terms of 
the distribution function, the particle stress energy ten- 
sor is (e.g. Ref. [36| ) 



t u » = / /pVdp, 



(187) 



expressed here in terms of the lab frame coordinate basis 
(note p^ — L^^p^). This stress energy obeys 



V M T^ = J p»C{f]dp. 



(188) 



In an astrophysical simulation involving neutrino trans- 
port, conservation of total energy-momentum is ex- 
pressed 



1 fluid 



E 

a=e,fj,,r 



T!; 



= 0, 



(189) 



where T^ id is the fluid stress-energy, and the sum is 
over all flavors of neutrinos. In light of Eqs. (|188l) and 
(|189p , the fluid energy and momentum source terms again 
emerge from momentum space integrals of collision inte- 
grals: 



- 7? V T v ^ — C ~ 
«iy v M J fluid — j^fj,q , 



where 



E 

a=e,fj,,r 



(190) 
(191) 

(192) 



(The energy and momentum equations are given by the 
projections of the stress-energy divergence orthogonal 
and tangent to the spacclikc slice, via contraction with 
the normal to the spacelike slice and the orthogo- 
nal projector 7^ = g^ u + n^n^.) Therefore knowledge 
of the distribution functions obtained by solution of the 
number-conservative Eq. (1178)) also yields the fluid en- 
ergy and momentum sources, via Eqs. (|190l) - (|192l) . 

Finally, we point out that while Eqs. (TTgg)) and (TIM)) 
are redundant with Eq. (|178l) and therefore do not need 
to independently be solved, they illuminate two chal- 
lenges to maintenance of global four-momentum conser- 
vation in numerical simulations. 

First, while the left-hand sides of Eqs. (Tl83l) and (fT84j) 
are transparently related to surface integrals in phase 
space, the presence of extra connection coefficient terms 
in the divergences appearing in Eq. (|189[) prevent trivial 



machine-precision four-momentum conservation in dis- 
cretized models. These extra source terms — Gt,k and 
{Gr,-y)j m Eqs. (|183[) and (|184[) . and analogous terms 
on the left-hand sides of Eqs. (I190[) and (|191[) — represent 
energy and momentum exchange between the forms of 
stress energy present (fluid and neutrinos) and the grav- 
itational field as embodied in the spacetime geometry. 
Even where global energy and momentum can be defined 
in general relativity — for instance in asymptotically flat 
spacetimes — conversion to a local conservation law, via 
definitions of gravitational energy and momentum den- 
sities and fluxes, is complicated at best (one possibility 
being use of a gravitational stress-energy pseudotensor 
0). 

Second, even if the gravitational source terms were 
not present in Eqs. (|183p and (|184[) . a straightforward 
discretization of Eq. (jl 78[) would not be consistent with 
straightforward discretizations of Eqs. (|183[) and (I184p . 
Presentation and testing of a full discretization is be- 
yond the scope of this paper. But as was done by 
Liebendorfer et al. [26[ in the spherically symmetric 
case, analysis of the relationship between the number- 
and four-momentum conservative formulations examined 
in Sec. lIIIEl can illuminate the search for discretizations 
that provide good energy conservation at modest res- 
olution in addition to the precise lepton number con- 
servation that comes naturally via Eq. (I178p . A few 
initial comments on this sort of consistent discretiza- 
tion of the spatially multidimensional angular moments 
formalism were also made in Appendix D of Ref. [35j . 
Consistent discretization of the spatially multidimen- 
sional Boltzmann equation will be algebraically more 
complicated than cither of these other two cases, as 
can be seen for instance by considering the manipula- 
tions between Eqs. (|156p and (|16ip . Our decomposi- 
tion in Eq. (|122p of the momentum space transforma- 
tion P l i into a portion Q l parallel to p l , and a por- 
tion U l i orthogonal to p l , is a helpful step. This al- 
lows us to introduce, via Eq. (|C24|) . the comoving- frame 

timelike components (Rn)° and (On)° into Eqs. (|123p 
and (|124p . The presence of all four components facili- 
tates more transparent analysis of the relationship be- 
tween the number- and energy-momentum-conservative 
formulations of the Boltzmann equation, by allowing use 
of Eqs. (|103p - (|106p to project out the different pieces 
of Eqs. (|123p and (|124p relevant to the energy- and 
momentum-conservative equations. However, it is not 
clear to us whether it will be possible in the discretized 
case to impose satisfaction of Eq. (I163p . which follows 
from Eq. (|162p . Use of this nontrivial identity is the fi- 
nal step we take in arriving at Eq. (|161l) . The extent to 
which a discrete analogue of Eq. (|163l) deviates from its 
vanishing analytic value — and what, if anything, might 
be done about it — remain to be seen. 

In spite of these remaining challenges to energy- 
momentum conservation, the number-conservative 3+1 
general relativistic Boltzmann equation we present in 
Eq. (|178p is a significant practical step forward towards 
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numerical implementation of multidimensional relativis- 
tic Boltzmann transport with coordinate basis spacetime 
position components and comoving frame curvilinear mo- 
mentum space coordinates. 



Appendix A: The 'extra' momentum derivative term 



In this appendix we compute the 'extra' momentum 
derivative term fT>% (P l iT l ^fip v p fl ), i.e. the fourth term 
on the left-hand side of Eq. (j2T))) . Three factors in this ex- 
pression are affected by the momentum derivative: P l ~ ll 
which is the only one in which momentum space connec- 
tion coefficients come into play; and the factors p v and 
j/ 1 . In the momentum derivatives, p° must be considered 
a function of p l through Eq. ([3]), so that 



dp 



Pi 



dp> 



dp 1 {—Pd) ' dp 1 



= S 3 t 



(Al) 



are the derivatives of the time and space momentum com- 
ponents respectively. 
First we show that 



ViP'i = 0. 

The result follows quickly with the help of Eq. (jT9 
dP 



■DjP 1 ; = 



dp 1 

dP\ 
dp 1 

dP\ 
dp 1 

dP\, 



n^p\ 



p\p\^\ ph 



pi pkpl^ 



dp k J 

$ dPh 



Qpk 



dp* J dp 
where we used the identity 



*>S=0. 







Qpi 

dp k 



■P>,. + P\- 



i dPh 



dp k 



(A2) 

l: 

(A3) 
(A4) 
(A5) 
(A6) 

(A7) 



as well as Eq. (fT7| . 

Next we turn to the first term in Eq. ([35]). Using 
Eq. ([All, 



- dp - 

= f 



1 



-V l ~ n-r>^ 4- V 1 -- 



(-PoV ° A 



(A8) 



(A9) 



where we have added and subtracted r°A-p AI to obtain 



the second line. 



Focus on the first term of Eq. (|A9[) . Using the lowered- 
index version 



dX ^ wVvV 



(A10) 



of the second geodesic equation, the first term of Eq. (|A9|) 
is 



(-p&Y 



/ dpy 
(-Pd) dX 

f_d^_ 

(~Pd) dX ' 



(All) 
(A12) 



On the right-hand side we use Eq. and the fact that 
Eq. ([3]) implies 



= d( PA /) =2p A d/, 



(A13) 



where here "d" represents any derivative operator, and 
find 



(-Pd) 



^ Of, POP" 



f dp 6 



(-P6) 



■Pd 



dX 



(A14) 



/ dp 



(-PbY 



(-Pd)- 



P^ (A15) 
• /'J",,, //'/'"• (A16) 



We rewrite the right-hand side using Eq. (|Alj) to obtain 



(-Pd) 



(-Pd) ®P % 
(-Pd) d P 



We massage the right-hand side further using the fact 
that p° is a momentum (three-) space scalar, and find 



/ 



-Pd) 



ofiPor 



(~Pa) 



(A19) 
. (A20) 



The sign reversal obtained in the last line will prove im- 
portant. 

Turning next to the second term of Eq. (|A9[) , we show 
that 



Invoking Eq. (jT2"j) . we have 

= (r% + l 



dL% 



dx^ 



(A21) 

(A22) 
(A23) 
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Because the coordinate basis connection coefficients are 
torsion- free (i.e. symmetric in the lower indices), 



Til/ T^U 



I— g dx^ 



dx^ 



(A24) 



(cf. Eg. (|3"Tj)'). But using the matrix identities (e.g. 
Ref. Hf) 



and 



Tr(lnA) = ln(detA) 
together with Eq. ©, we also have 

dL v dln(-g)- 1 / 2 d\n^ 



dx^ 



dx^ 



(A25) 



(A26) 



(A27) 



and Eq. (IA21[) is confirmed. 

Combining Eqs. (|A20j) and (|A21|) in Eq. (fAUj) . we have 



dp 
dp 1 1 



fT- 1 — w# 

J L VfJ, Q _{ P 



-/rv A pVn(-Po)^ 



(-po) 



(A28) 



for the first term in Eq. (f35|) . 

Next we turn to the second term in Eq. (|35|) . which is 
more straightforward. Using Eq. (jAll) . 



r ^ dp* = 
= / 



1 



(-P0) 
1 



(A29) 



(A30) 



where we have added and subtracted r -^" to obtain 
the second line. The first term in Eq. (|A30[) vanishes, 
because 



(A31) 



with both terms in the middle independently vanishing 
thanks to Eq. (|A13|) . As for the second term in Eq. (|A30I) . 
we have 



L^uL^I/a r%„ + V^y^^JL (A32) 



/Li Z-' J - / /J, L Vp 



dxP 



dx» ■ 

Combining these results, we have 



/ r "MP apt -/f ^ ^+ dx » 

for the second term in Eq. (|35l) . 



(A33) 



(A34) 



Appendix B: Tetrad solution 

In this appendix we give the equations that com- 
plete the specification of the comoving frame tetrad 
w+ u+, ?/+ of orthonormal basis vectors, according to 
the algebraic conditions outlined in Sec. IIII Bl The time- 
like basis vector is the Lagrangian observer four- 
velocity, which is also the fluid four-velocity given by 
the evolution of the fluid and geometry. In Eqs. (J7DJ)- 
(|72|). the spatial basis vectors -u+,y M ,z M are decom- 
posed into portions orthogonal — An^ .Bn^ 1 : Cn^ — and 
tangent — a Ai ,6' i ,c At — to the spacelike slice. The orthogo- 
nal components A,B,C are given by Eqs. (I75)) - ([5D|) ; and 
in the lab frame coordinate basis, a M , 6 M , c M are spacelike 
(Eq. I|74p). so that the nine components a l ,b l ,c l remain 
to be determined. For starters, we specify 



b 1 = 0, 

c 1 = c 2 = 0. 



The single component c 3 is given by 

733 (c 3 ) 2 = 1+ (v 3 c 3 Y 



(Bl) 
(B2) 

(B3) 



which follows from Eqs. dHDJ), ([55]1. and (|B2]l . Once c 3 
is known, the non- vanishing components b 2 and b 3 are 
determined by Eqs. and (|90p . which become 

722 (b 2 ) 2 + 2 723 (b 2 b 3 ) + 733 (6 3 ) 2 = 1 + (v 2 b 2 + v 3 b 3 ) 2 

(B4) 

and 

723^ 2 c 3 + 733& 3 c 3 = (v 2 b 2 + v 3 b 3 ) (v 3 c 3 ) (B5) 

when Eqs. (j75 ]> - (|5D )1 . ([BT|) - (|B2|l are taken into account. 
Now that b l and c 1 are specified, the three components 
of a 1 are given by Eq. 



7n (a 1 ) 2 + 2 7i2 (Va 2 ) + 722 (a 2 ) 2 
+ 2723 (a 2 a 3 ) +733 (a 3 ) 2 

= 1 + (via 1 +v 2 a 2 + « 3 a 3 ) 2 , (B6) 



by Eq. 

12 2 2 3 2 

712a b + 722a b + 7 32 a b 
+7i 3 a 1 6 3 + 7 23 a 2 6 3 + 733a 3 & 3 

= (via 1 + u 2 a 2 + v 3 a 3 ) (v 2 b 2 + i^b 3 ) , (B7) 

and by Eq. 



7i 3 a 1 c 3 + 723a 2 c 3 + 733a 3 c 3 



= (via 1 + v 2 a 2 + v 3 a 3 ) (v 3 c 3 ) , (B8) 
where Eqs. ([7g ]) -([8"0] ) . ([BT] ) -([E2] ) are taken into account. 



Appendix C: Comoving frame connection 
coefficients 

In this appendix we derive, in the context of the 3+1 
approach, expressions for the comoving frame connection 
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coefficients Y % op, — and their contraction with the momen- 
tum variable transformation P l j — appearing in the mo- 
mentum space divergence in Eq. (j42|) . We begin with the 
expression 



VP - ~ = L? LP - fV L"M 



(CI) 



from Eq. (|118p . Key to our derivation is use of the Eule- 
rian decompositions of the tetrad L M A and its inverse LP ^ 
given in Eqs. (p?3"|) and keeping in mind that their 

projections £p a and ^ are tangent to the spacelike slice 
(see Eqs. (J95J), CUT]), ((Ml), and ®). 

Some relations involving derivatives of the unit normal 
n,P will prove useful. The gradient of the unit normal is 
related to the extrinsic curvature and lapse function by 
401 



rifx da 
a dx v 



Because 



n v V^n v = V» {n v n v )/2 = 0, 



(C2) 



(C3) 



and because K^ v is tangent to the spacelike slice, i.e. 
spacelike in the lab frame coordinate basis, 



nPR,, 



n v K„ 



0. 



the nonvanishing projections of Eq. (|C2[) are 

, m 1 da 

— — — , 

a ox v 



(C4) 

(C5) 
(C6) 



in the lab frame coordinate basis. Equation (|C5|) relates 
the four-acceleration of an Eulerian observer to the gra- 
dient of the lapse function. Equation (|C6|) relates the 
spatial part of the gradient of the unit normal to the ex- 
trinsic curvature, expressing the fact that the direction 
of the normal varies with the warp of the slice as embed- 
ded in spacetime. Another relation valid in the lab frame 
coordinate basis for vectors zP tangent to the spacelike 
slice (z° — 0) is 

This follows from writing 



dn.P dn° dri 1 



(C8) 



= *»*B? + «B? (C9) 



and using goi — ft and Eq. (J51 

We proceed with the calculation in three parts, based 
on the terms 

V^i"* = Co (V M n w ) + n v {d^C ) + (V^f e ) (CIO) 

obtained by using the Eulerian decomposition in Eq. (|93[) 
in the rightmost factor of Eq. (IC1[) . 



We consider first the contribution 

LP U LP^ Co (V^n") = IK (n"£ A + Pp) C (V^n v ) 

(Cll) 

to Eq. (|C1[) from the first term in Eq. (IC10I) . where we 
have used Eq. f|C3j) . Using Eq. (|C2j) we have 



for this contribution. 
Next is the contribution 



(C12) 



LKL^P, n v (d^Co) = -C p (n%> + *" A ) (d^Co) (C13) 



to Eq. (jCip from the second term in Eq. (IC10I) , where we 
have used the fact that i p v is spacelike. Using Eq. (|55p . 
we have 



L p v L^n v (d^Co) 
_CP_ 
a 



Of 



dx % 



(C14) 



for this contribution. 
Finally, the contribution 



ifivWp. (V M r e ) = {C p n v + (n"£ A + *" A ) (V^ fi ) 

(C15) 

to Eq. (IC1[) from the second term in Eq. (jClOl) is more 
complicated. Using the identity = V M in v i v o) — 
(V M n„) l v 6 + n v {V^to) along with Eq. jC2j, two of the 
terms are 

(C p n v ) (nP £ A + -P A )(V M f») 

= -^f^(r l %+f A )(V (1 n 1/ ) (C16) 



Another term gives 



-l pv CfitpD {V v n p -d v n p ) 



(C18) 
(C19) 

(C20) 



where now Eq. (|C7[) has been used in addition to 
Eq. ([C2]l . The final piece is 



t pk p , 



( d£ kC , _ £ l o djj. 
V dxi 2 dx k 



(C21) 
(C22) 
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in which we have anticipated that the two terms from the 
connection coefficients antisymmetric in i and j will can- 
cel upon contraction with the product (? (which be- 
comes symmetric in i and j upon contraction with p^p"). 

In summary, the comoving frame connection coeffi- 
cients TPpf, are giv en by the sum of Eqs. (|C12j) . (|C14j) . 
(|CT7)) . (fU2T)|) . and (IC221) . We note that two of the terms 
involving Kjk cancel upon contraction with p^p v . 

In the momentum space divergence in Eq. (|42l) . the 
comoving frame connection coefficients appear con- 
tracted with the momentum space coordinate transfor- 
mation P % % and two factors of momentum. It will prove 
convenient to use the decomposition of P l i in Eq. (If 07[) : 



which follows from Eqs. © and ([JT3]) . Thus Eq. (jC23|) 
becomes 



PhT^ap^ 



Q 1 ( Pq) v 6 , 771 r J 



(C25) 

It is in this combination that we use the comoving frame 
connection coefficients — the sum of Eqs. (ICI2I) . (ICI4I) . 
(|CI7|) . (|C20|) . a nd (|C22j) — in the momentum space di- 
vergence in Sec. IIII D[ 



p\r iiil p v p» = 



Pi 



p 



+ U l l )r 0ii p lJ p^ (C23) 
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